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ABSTRACT 

We present high resolution H-band polarized intensity (PI; FWHM = O/'l: 
14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar 
disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance 
of 28 AU (0."2) up to 210 AU (l."5). In both images, a giant inner gap is clearly 
resolved for the first time, and the radius of the gap is ~70 AU. Our data show 
that the geometric center of the disk shifts by ~6 AU toward the minor axis. 
We confirm that the brown dwarf companion candidate to the north of PDS 70 
is a background star based on its proper motion. As a result of SED fitting by 
Monte Carlo radiative transfer modeling, we infer the existence of an optically 
thick inner disk at a few AU. Combining our observations and modeling, we 
classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the 
analysis of L'-band imaging data, we put an upper limit mass of companions at 
~30 to ~50Mj within the gap. Taking account of the presence of the large and 
sharp gap, we suggest that the gap could be formed by dynamical interactions 
of sub-stellar companions or multiple unseen giant planets in the gap. 

Subject headings: planetary systems — protoplanetary disks — stars: individual 
(PDS 70) — stars: pre-main sequence — polarization 



1. Introduction 



Protoplanetary disks are believed to be the birthplaces of planets (e.g.. lHayashi. Nakazawa. &: Nakagaw 



19851 1: hence, understanding the evolution of these disks guides our understanding of the pro- 



cess of planet formation. Disks which have substantial infrared excesses but reduced fluxes 
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27 Kavli Institute for the Physics and Mathematics of the Universe, The University of Tokyo, Kashiwa 
227-8568, Japan 

28 Dcpartment of Cosmosciences, Hokkaido University, Sapporo 060-0810, Japan 
29 Astronomical Institute, Tohoku University, Aoba, Sendai 980-8578, Japan 

1 Based on data collected at the Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 
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at wavelengths <20 /im, i.e., transitional disks (jStrom et al.lll989f). could be relat ed to the 



early phases of planet formation (see a recent review of I Williams fc Ciezal 1201 if ) and are 



therefore particularly important for understanding how, where, and when planets form. For 



terferometry at (sub) millimeter wavelengths 



(e.g., 


Pietu et a 


. |2006 




Andrews et al. 


2011) 


(Fukagawa et al. 




20061 


Thalmann et al. 


2010 



Hashimoto et al.l 1201 ll ) . Numerous mechanisms have been proposed to explain the clear 
ing of gaps in transitio nal disks, including grain growth (e.g ., iDullemond fc Dominikll2005l ) 



Papaloizou et al 



2007 



Clarke. Gendrin. fc Sotomavorll200l!). and g ravitational interactions 



Zhu et al.l 120 111 ). Two possible methods 



photoevaporation (e.g. 
with orbiting planets (e.g. 
to distinguish the disk-planet interactions from other aforementioned proposed gap-clearing 
mech anisms could be the det ection of (1) a planetary companion in the inner hole/gap region 
(e.g., iKraus fc Ireland! |2 0121 ) or (2) a ring-like gap b etween optically thick inner and outer 
disks (i.e., pre-transitional disk; lEspaillat et al.l 120071 ) because dynamical formation of wide 
gaps could be the on ly surviving mechanism for wide gapped disks (e.g., iPapaloizou et al. 



20071 : IZhu et al.ll201lh 



One good candidate to investigate the inner hole/gap region at tens AU in the disk in 
pre-transitional disks is the weak-lined T Tauri sta r PDS 70 (K5 type; 0.82M^; <10 Myr ; 
Gregorio-Hetem fc Heteml I2OO2I ; llliaud et all bood : iMetchev. Hillenbrand, fc Meverl I2004T ). 
A scattered light disk with a radius at 14 to 140 AU was detected by Ks-baiad imaging 
(IRiaud et al.ll2006l ). The possible presence of inner and outer di sks w ith different t emper - 
atures were suggested by IMetchev. Hillenbrand, fc Meyerl (120041 ) and IRiaud et al.l (120061 ). 
which may imply that PDS 70 is a pre-transitional disk object. In this Letter, we present 
high resolution imaging of PDS 70 with Subaru/HiCIAO and Gemini/NICI. 



Observations & Data Reduction 



2.1. iJ-band Polarimetry with Subaru/HiCIAO 



Observations of PDS 70 were conducted with HiCIAO (ITamura et al.l 120061 ) on the 
Subaru telescope in polarized d ifferential imaging (PDI) mode, combined with angular dif- 
ferential imaging (ADI) mode (IMarois et al.ll2006l ). on 2012 February 28 UT. H-band lin- 
ear polarization images were taken under the prog ram SEEDS (Strategic Explorations of 
Exoplanets and Disks with Subaru; ITamura 1 120091 ). In PDI+ADI, we employed a double 
wollaston prism to split incident light into four images, each has 5" by 5" field of view with 
a pixel scale of 9.5 mas/pixel, to make the saturated radius as small as possible. 
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Polarization was measured by rotating the half waveplate to four angular positions (in 
the order of 0°, 45°, 22.5°, and 67.5°). We obtained 22 full waveplate rotation cycles, taking 
a 15-s exposure per waveplate position. The total integration time of the PI image was 
660 s. Even wit h a large airmass of ~2 during observations of PDS 70, the adaptive optics 



system (A0188; lHavano et al.ll2004h provided a stellar PSF of PDS 70 with FWHM of 0."1 



in the H band. The angle of the total field rotation was ~13°. 



The polarimetric data were reduced with the same procedure as for lHashimoto et al. 



(boilh using IRAfB Since the Strehl ratio of the stellar PSF of PDS 70 was ~0.26, a stellar 



halo remained in the image, i.e, the PSF convolved by seeing was not perfectly corrected by 
the AO and partially remains. This halo appears to be polarized along the minor axis of 
disks since the polarization of forward and back scattering is smaller due to the deviation 
from 90° in the scattering angle, and then the stellar halo has a significant net-polarization. 
In other words, since the flux of a PSF convolved by seeing contains contributions from both 
the central star and disk, the seeing-PSF might be polarized along the minor axis of a disk. 
To correct for this "polarized halo", we first derived the net-polarizarion of PDS 70 with 
aperture polarimetry. We then made a model "polarized halo" of PDS 70 by multiplying 
the added image of o- and e-rays by the derived polarization of P = 0.503% ± 0.001% with 
9 = 66.32° ± 0.06°. The errors were calculated with the photometric errors, and are thus a 
lower limit. Finally, we subtracted the model "polarized halo" and obtained final Stokes Q 
and U images. Figure [1] demonstrates the subtraction of the "polarized halo" . 



2.2. L'-band imaging with Gemini/NICI 

We carried out ADI observations of PDS 70 in the V band on 31 March 2012, using 



the N ear-Infrared Coronagraphic Imager (NIC!) and the 85-element AO system (IChun et al. 



20081 ) mounted on the 8.1 m Gemini South telescope at Cerro Pachon, Chili. NICI utilizes 
a 1024 x 1024 ALADDIN InSb array with an plate scale of 18 mas/pix and field of view of 
18 x 18 arcsec. 

We obtained 145 frames for PDS 70 with 0.76-s x 25 coadd for V band. Although the 
sky was clear for the full night, the seeing was variable, therefore a subset of 95 frames were 
combined to produce the final image. The total exposure time was 1805 sec and the spatial 
resolution was achieved 0."11. The angle of the total field rotation is ~100°. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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The data reduction was performed with the LOCI algorithm ( Marois et al.ll2006l ; lLafreniere et al. 



20071 ) in annular regions of 1000 x FWHM (JVa = 1000) with reference images selected from 
frames with at least 0.5 x FWHM field rotation (Ng = 0.5 ). A large optimization area 
dimin ishes the impact of the disk on the optimization process fjBuenzli. Thalmann. fc Vigan 



2010l ). The ratio of radial and azimuthal width of the optimization area (g) is 1. 



Results & Discussion 



3.1. H-band Polarimetry and L'-band imaging 

Figure [2] shows LL-band PI images and the L'-band LOCI image of PDS 70 assuming a 
distance of 140 pc, along with radial surface brightness profiles. We find a clear e lliptical ring 



in the if-band, which has not been reported in previous high-resolution imaging fjRiaud et al. 



20061 ). A partial elliptical disk is observed in the L'-band, due to the inevitable loss of flux 
in the process of LOCI; hence, we derive radial profiles of suface brightness based on the 
LLband PI image only and companion mass limits from the L'-band LOCI image only. 

We consider that the ellipse shape is due to the system's inclination, and show the 
results of fitting an ellipse to these data in Table [TJ The position angle of the major axis 
and the inclination of the disk are similar with those of Lf s -band imaging (PA w 155° and 
i « 62°; iRiaud et aliboOfih . We measured an offset of 44 ± 3 mas (~6 AU) at PA=87.9° 
between the geometric center of the disk and the central star. The positional accuracy of 
the central star is 1.5 mas (0.2 AU). The direction of this offset is roughly consistent with 
that of the minor axis of 68.6°, and the sign of this offset indicat es that the southwe st side 
is inclined toward us (i.e. the near side) (see the model image in iDong et al.ll2012bl ). This 
geometry is also consistent with the f acts that (1) the northe ast side of the disk is wider due 
to the back illumination of the wall (IThalmann et al.l 1201011. and (2) the southw est side is 
brighter than the northwest side due to forward scattering (IFukagawa et al.ll2006l ). 



Assuming the cavity edges correspond to the peak PI of the disk, a radius of the cavity is 
measured as ~70 AU. The outer radius of the disk is measured to be ~140 AU in figure [2(d) , 
and corresponds to the location at which our sensitivity is no longer sufficient to detect 
extended emission. 

A single power-law fit was performed to the radial profiles along the minor and ma- 
jor axes (figure EH and e). Our results of ~ -3.6 and ~ -1.7 are different from ~ -2.8 of 
Riaud et all hood ). 



We also found a flux deficit of PI in the direction of the minor axis. Since the observed 
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scattering angle at the minor axis deviates from 90°, the polarization fraction along the 
minor axis is lower. The PI at the minor axis is therefore lower than at the m ajor axis, and 
such hole-like structures are similar to those discussed in iPerrin et al.l (120091 ) . 



We c hecked the proper m otion of the companion cadidate to the north of PDS 70 re- 
ported by iRiaud et al.l (I2006|). PDS 70 has a proper motion of ( uncos 5, fj,$) = (—24.7 ± 
11.4, — 13.3 ± 11.4) mas/yr ( iRoeser. Demleitner. fc Schilbachl 120101 ); the separation between 
PDS 70 and the companion cadidate should increase if the companion cadidate is a back- 
ground star. Since the separation in HiCIAO a nd NICI images ar e 324.13 ± 0.15 AU and 
324.44 ± 0.10 AU, respectively, the separation in lRiaud et al.l ( 120061 ) is inferred to be 309±12 
AU. Our e stimation has a goo d agreement with the actual observed separation of 301.75 ± 
0.06 AU in IRiaud et al.l ( 120061 ). and therefore, we concluded that the companion cadidate is 
a background star. 



3.2. Detectable planetary-mass companions 



Since the follow-up L'-band observations with Gemini/NICI failed to detect any signif- 
icant signals of point-like sources in the gap we put constraints of upper limits for compan- 
ion^). Figure |2Jf) shows the detectable masses of companions at 5cr. The LOCI parameter 
of the optimization area is 250 x FWHM, which is different from that described in sec. 12.21 
For that, we first applied a median filter with . " 1 1 width to the image, and then calculated 
the standard deviation as a noise level in concentric annuli along the m ajor axis of the disk . 
The mass was calculated by assuming the COND evolutionary model (iBaraffe et al.ll2003l ). 



a distance of 140 pc, and an age of 10 Myr (IMetchev. Hillenbrand, fc Meyerll2004[ ). We took 



into account the flux loss due to the partial self-subtraction by testing how point sources are 
affected by LOCI. The detectable mass limit is tens Mj, therefore, stellar companions down 
to masses associated with massive brown dwarfs are excluded within the gap. 



3.3. Modeling of spectral energy distribution (SED) 



Although the SEP fitting for PDS 70 has been performed in previous studies ( IMetchev. Hillenbrand, fc 



20041 ; IRiaud et al.ll2006l ). both the availability of new archival photometric data (table [2]) and 



our imaging results motivates us to revisit the SED of the system using Monte Carlo radia- 
tive transfer (MCRT). Note that though PDS 70 has been observed by Spitzer IRS, the 
object was mispointed by ~2.3" with components along both the spatial and dispersion 
axes, therefore, the IRS spectrum of PDS 70 is not used in this work. 
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Se tup for modeling;. The m ethod for our MCRT simulations is desc ribed inlDong et a 
f l2012a| ) and IWhitney et al.l (120121 in prep.). In a subsequent paper II flDong et al.ll2012b 
in prep.), we will perform a detailed radiative transfer modeling of both the SED and the 
SEEDS imagery, and present a fiducial disk+cavity model which reproduces both observa- 
tions well. We will also explore the parameter space around the fiducial model in paper 
II, and provide a full discussion on the constraints on the various model parameters there. 
In this letter, we only briefly describe the fiducial model and its resulting SED. We note 
that this model is not a best fitting model in an absolute sense, since a full x 2 fitting of the 
observations with all the free param eters in a protoplanetary disk is essentially impossible 
f lMathews. Williams fc Menardll2012[ ). However the constraints on many parameters such as 
the cavity size, depletion, and the surface density of the inner disk are reasonably tight, as 
will be shown in paper II. 

Our model contains a cavity 70 AU in radius. The surface density both inside and 
outside the cavity decreases with radius as £ oc ^-e~ R ^ Rc , where R c = 50 AU, while E inside 
the cavity is reduced to 5x the extrapolated value from the outer disk, with 5 being the 
depletion factor. The temperature structure of the disk is determined from the radiative 
transfer calculations. The inner edge of the disk is self-consistently determined at the dust 
sublimation temperature (~ 1600 K). We ignore accretion in the model , as suggested by its 
nature of being a weak line T Tauri star (jGregorio-Hetem fc Hetemll2002l ). We use a pre-main 
sequence star of spectral type K5, radius 1.39 R^, mass 0.82 M ffl an d te mperature 4500 K fo r 
the central source, as suggested by lGregorio-Hetem &: Heteml (120021 ) and lRiaud et al.1 (120061 ). 
The disk has a gaussian density distribution in the vertical direction, with scale heights h 
as input parameters. Two disk components are included: a thick disk with small grains 
(sub-micron-sized), representing the pristine ISM-like dust; and a thin disk with large grains 
(up to ~mm-sized) and 20% of the scale height o f t he sm all grains, as the result of grain 
growth and settling ( Dullemond fc Dominik 2004al lbl. 12005 ). Most of the dust mass (0.967) 
is in the settled disk, and the total dust-to-gas ratio is assumed to be 1/100. The SED is 
produced assuming a disk inclination angle 50°. The other parameters are summarized in 
table m 



SED fitting. Figure [3] shows the good agreement between our model SED and available 
photometric data. As we will show in paper II, the thermal emission from the cavity wall at 
~70 AU p eaks at ~ 40 /xm. The wall emission is a major signature of (pre-)transitional disk 
SED 



e.g. 



Espaillat et al.ll2007f ). 



The surface density of small dust is ~ O.OOlg cm -2 at 0.1 AU, and the opacity of small 
dust is ~ 10 4 cm 2 g" 1 at ~ 1/im (roughly the peak of the stellar radiation). This makes the 
inner disk optically thick in the vertical direction. Since we assume a surface density profile 
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decreasing with increasing radius, when moving out the cavity gradually becomes vertically 
optically thin, and the transition happens at ~ 1 AU if only taking into account the small 
dust (the existence of big dust inside the cavity is poorly determined from SED and scattered 
light image). As we will discuss in detail in paper II, this optically thick inner disk is needed 
to explain the ~2-40 /im SED, as models with op tically thin inner disk fail to reproduce the 
SED. Following the convention in the literature (jEspaillat et al.l 120071 ). we classify PDS 70 
as a pre-transitional disk object. 



3.4. Origin of the gap 



Grain growth is one mechanism which can potentially form disk gaps (e.g.. lBirnstiel. Andrews, fc Ercol< 



20121 ). It is capable of reproducing the SED of transitional disks but not observed millimeter- 
wavelength images. However, the sharpness of PDS 70's if-band gap edge suggests grain 
growth is unlikely to be the primary reason, since the m odels generally predict smooth gap 
edge features (e.g.. iBirnstiel. Andrews. &: Ercolanoll2012l ). 



Photoevaporation is another possible mechanism. A photoevaporative wind can prevent 

outer disk material from feeding the inner disk; without this suppl y, the inner disk material 

will r apidly accrete onto the central star, creating a cavity inside out (IClarke. Gendrin. fc Sotomayor 
20011 ). Although a pre-transitional disk-like structure c ould be produced for a very short pe- 
riod (see figure 1 in lAlexander. Clarke, fc Pringldl2006l ). it is quite unlikely to observe such a 
snap-shot by coincidence, so that it may be difficult to explain the existence of the optically 
thick inner disk inferred from the SED. 



Dynamical interacti ons with (sub)stellar companions (e.g.. lArtvmowicz fc Lubowlll994l ) 



or orbiting planets (e.g., iPapaloizou et al.ll2007l ; IZhu et al.l 120111) are also po tential mecha 



nisms, as several binary systems (e.g., CoKu Tau 4; llreland fc Krausl 120081 ) are known to 
have a cavi ty in their circum binary disk, and simulations show that multiple planets can open 
wide gaps (IZhu et al.ll201 11). Furthermore, the survival of an optically thick inner disk within 
a few AU in lZhu et al.l (120 llh mimics the inner disk of PDS 70. Our L'-band observation has 
only ruled out companions with a mass of >50Mj in the gap region. Therefore, dynamical 
interactions with low mass brown dawarfs or giant planets may be the origin of the gap. To 
put a futher constraint of an upper limit mass, future obs ervations should be pu rsued, such 
as the aperture-masking interferometric observations (e.g.. iNakajima et al.lll989l ). which has 
a better contrast than LOCI imaging. Since PDS 70 harbors the pre-transitional disk with 
one of the largest gaps, this object may be one of the best candidate for future high-contrast 
planet imagers. 
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Table 1: The results of the ellipse and the SED fitting for the disk of PDS 70 



Ellipse fitting 1 



Diameter of the major axis 2 (AU) 


137.0 ± 1.1 


Diameter of the minor axis 2 (AU) 


88.6 ± 0.5 


Position angle of the major axis (°) 


158.6 ± 0.7 


Inclination by ellipse fitt 


ing 3 (°) 


49.7 ± 0.3 


Geometric center 2,, (AU : 


AU ) 


(6.1 ± 0.3 : 0.2 ± 0.4) 


SED fitting 5 


M disk 




0.003 


/ 




0.967 


/i°(100AU) 




2 


h° 
°b 




1.2 


/i£(100AU) 




2 


K 




1.2 






10" 3 


/i°(100AU) 




10 






1.2 


^(lOOAU) 




10 


K 




1.2 


"cav,s 




10~ 3 



1 In the ellipse fitting for the disk, the peak positions were first directly determined by the radial profile at 
position angles every 5° the disk. We then conducted an ellipse fit using an implementation of the nonlinear 
least-squares Marquardt-Lcvenbcrg algorithm with five free parameters of lengths for the major and minor 
axes, position angle, and central positions. 

2 Wc assume that the distance of PDS 70 is 140 pc. 

3 Derivcd from the ratio of the major and minor axes. 

4 Central position (0, 0) is corresponding to the stellar position. 

5 Row (1): Total mass of the disk (assuming a dust-to-gas-ratio 100). Row (2): Mass fraction of big dust in 
total dust. Rows (3), (5), (8), and (10): Scale height h at 100 AU. Subscripts b and s indicate big dust and 
small dust, respectively. Superscripts o and c indicate outer disk and cavity, respectively. Rows (4), (6), (9), 
and (11): Power index b in h oc R b in various disk components. Row (7) and (12): Depletion factor of the 
big and small dust disk. 



Table 2: Archival photometry data for PDS 70 



Wavelength 


F„ (mJy) 


Note 


Plotted color 
in fig. [3] 


ua ,b 


9.7 


Greeorio-Hetem et al. (1992) 


cyan 




41.8 


Greeorio-Hetem et al. (1992) 


cyan 


ya,b 


99.2 


Greeorio-Hetem et al. (1992) 


cyan 


R a ' b 


160.9 


Greeorio-Hetem et al. (1992) 


cyan 


ja,,b 


216.3 


Greeorio-Hetem et al. (1992) 


cyan 


2MASS (J) a < b 


311.9 ± 6.9 


Cutri et al. (2003) 


magenta 


2MASS (HY' b 


342.7 ± 12.6 


Cutri et al. (2003) 


magenta 


2MASS (Ks) a ' b 


275.6 ± 5.8 


Cutri et al. (2003) 


magenta 


WISE (3.4 /mi) b 


188.1 ± 4.0 


Cutri et al. (2012) 


green 


WISE (4.6 /mi) b 


142.0 ± 2.6 


Cutri et al. (2012) 


green 


WISE (12 /im) b 


153.9 ± 2.3 


Cutri et al. (2012) 


green 


WISE (22 /im) b 


341.8 ± 0.7 


Cutri et al. (2012) 


green 


AKARI (9 fjon) 


201.2 ± 25.8 


VizieR 11/297 


black 


AKARI (18 £tm) 


209.8 ± 13.4 


VizieR 11/297 


black 


AKARI (90 /tm) 


851.1 ± 62.6 


VizieR 11/298 


black 


IRAS (12 //m) 


251 ± 25.1 


Moshir (1989) 


blue 


IRAS (25 /iin) 


348 ± 27.8 


Moshir (1989) 


blue 


IRAS (60 /iin) 


884 ± 61.9 


Moshir (1989) 


blue 


MIPS (24 /iin) 


349.7 C ± 7.0 


Spitzer Heritage Archive 


red 


MIPS (70 /iin) 


1049.9 C ± 19.9 


Spitzer Heritage Archive 


red 


MIPS (160 /iin) 


873.0 c ' d ± 36.2 


Spitzer Heritage Archive 


red 



'An extinction law was adopted from Mathis (1990). 
' Absolute flux conversions i n optical, 2MASS, 



WISE pho tomet ric data w ere adopted from 



Bessell. Castelli. k Pled (|1998l ). ICohen. Wheaton. fc Megeathl (|2003h . and Uarrett et all (|201lh . respectively. 



: Sincc photometric data were not available, we conducted aperture photometry for archival images. 

'Since the sky value was not sufficiently measured due to a small field-of-view of available data, our photometry 

may be less reliable. 

'This work is based in part on observations made with the Spitzer Space Telescope, obtained from the NASA/ 
IPAC Infrared Science Archive, both of which are operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with the National Aeronautics and Space Administration. 



„ v / \ 

(b)/-- 




Fig. 1. — i^-band polarization vectors of PDS 70 are superposed on the PI image with a 
software mask with 0."4 diameter before subtracting polarized halo (a) and after subtraction 
(b). The plotted vectors are binned with spatial resolution. The field of view (FOV) is 3."0 
x 3."0. All plotted vectors' lengths are arbitrary for the presentation purpose. 
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(d) //-band PI Surface brightness (e) //-band PI Surface brightness (f) L '-band LOCI companion limits 




Radius (arcsecond) 



Fig. 2. — (a): fZ-band PI image of PDS 70 with a software mask with 0."4 diameter, (b): 
same with (a), but its features. The solid ellipse indicate the ring- like disk. The filled circle 
represent the geometric center of the disk, (c): L'-band LOCI image of PDS 70 with a soft- 
ware mask with 0."4 diameter. The parameters in LOCI reductions are described in sec. 12.21 
The FOV of three images are 3."0 x 3."0 with a convolution of a spatial resolution, (d) and 

(e) : Radial profiles at yellow hatched regions of minor and major axes in (b). The values of 
the profile at northwest and southwest are multiplied by ten for the presentation purpose. 

(f) : Detectable mass at 5<r based on the L'-band LOCI image. The LOCI parameters are 
same with those described in sec. I2.2[ but the optimization area is 250 x FWHM (A/a = 
250). 
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Wavelength [jim] 

Fig. 3. — Pre-transitional disk model of PDS 70. Filled circles represent archival photometry 
(see table [2] for photometry data). The solid black line is the best-fit model with a gap of 
~70 AU (see table Q] for model parameters). Separate model components are as follows: 
stellar photosphere (green dotted line), thermal emission (blue dotted line), and scattered 
light (red dotted line). 
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